function y = computeGc(element,B)
PosDomain = element.PosDomain();
if strcmp(element.mode,'constant')
    y = zeros(3,2);
    for i = 1:size(PosDomain(:),1)
        shift = 2*(PosDomain(i)-1);
        y(1,1) = y(1,1) - B(1,1 + shift);
        y(2,2) = y(2,2) - B(2,2 + shift);
        y(3,1) = y(3,1) - B(3,1 + shift);
        y(3,2) = y(3,2) - B(3,2 + shift);
    end
    J = computeJ(element);
    y = y*J;
elseif strcmp(element.mode,'linear')
    % TODO
end

end
